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We elaborate further on the possibility that the inflationary primordial power spectrum contains 
superimposed oscillations. We study various effects which could influence the calculation of the mul- 
tipole moments in this case. We also present the theoretical predictions for two other cosmological 
observables, the matter power spectrum and the EE polarization channel. 
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INTRODUCTION 



II. ROBUSTNESS OF THE WIGGLES 



The possibility that the multipole moments Ci char- 
acterizing the angular distribution of the Cosmic Mi- 
crowave Background Radiation (CMB) anisotropy on the 
celestial sphere possess superimposed oscillations has re- 
cently been investigated in Ref. pj. In that article, the 
superimposed oscillations originate from trans-Planckian 
wiggles [2] in the primordial inflationary power spec- 
trum but the study of Ref. [l| was meant to be as in- 
dependent as possible from the details of the underlying 
model. The presence of oscillations in the primordial 
spectrum has also been envisaged in Refs. jjj. Then, 
it has been demonstrated that the superimposed oscil- 
lations can cause a significant improvement of the fit 
to the first year Wilkinson Microwave Anisotropy Probe 
(WMAP) data ||, thanks to the presence of the cosmic 
variance outliers around the first Doppler peak. More- 
over, it has also been shown that the corresponding drop 
in the x 2 is statistically significant according to the so 
called F test. We have since received various inquiries 
about different effects that could modify the result ob- 
tained in Ref. Q] like the influence of the splinning, the 
consideration of the lensing [f| and the way of estimat- 
ing the statistical significance of the drop in the \ 2 ■ In 
this addendum, we clarify these issues and, in addition, 
present new results which are important for the com- 
pleteness of Ref. namely a fit with an improved x 2 
which does not suffer from the back reaction problem and 
the prediction for two other cosmological observables, the 
matter power spectrum and the EE polarization chan- 
nel. Finally, we would like to emphasize that the general 
questions analyzed in this Brief Report are important ir- 
respective of the data set used to study them. Therefore, 
although we use the first year WMAP data, the con- 
clusions reached in this addendum are also valid for the 
future CMB data releases. 
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In this section, we investigate effects that could possi- 
bly influence cither the numerical calculation of the CMB 
multipole moments or the statistical analysis of the sig- 
nificance of the "wiggles" . 



A. Multipoles splinning 

In the case where oscillations are present in the ini- 
tial power spectrum, a correct numerical computation of 
the CMB multipole moments requires one to significantly 
boost the numerical accuracy of the code used to derive 
the so-called transfer functions |(|. Indeed, it is necessary 
to know these functions with a high precision if one wants 
to correctly transfer the contribution of the primordial 
oscillating spectrum into each multipole moments. In 
Ref. 0], such high accuracy computations have been per- 
formed with the CAMB code [7| and various tests have been 
carried out in order to ensure that the multipole moments 
were properly computed. However, one must also pay at- 
tention to the sampling and the splinning performed on 
the multipole moments. They are both present by default 
to avoid prohibitive computation time. The sampling re- 
quires the computation of some multipole moments only 
and is scale dependent (i.e. ^-dependent) while the splin- 
ning interpolates between those multipole moments and 
thus allows to recover the complete angular spectrum. 
In Ref. |IJ , the sampling and splinning have been kept to 
their default option which is clearly not very appropriate 
to the case where superimposed oscillations are present. 
Indeed, as long as the multipole moments Ci are not 
computed for each £, there is the danger to undersample 
the signal which could result in the appearance of oscil- 
lations with an incorrect frequency. At large scales, the 
multipole moments are always calculated for each value 
of I and at very small scales, the sampling has no effect 
since the oscillations considered in Ref. [IJ are logarith- 
mic in the Fourier space. Therefore, the danger is partic- 
ularly present at intermediate scales, i.e. around the first 
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Doppler peak. In this Brief Report, we have removed the 
sampling and the splinning from the CAMB code such that 
each multipole is now computed. One finds that this does 
not affect the determination of the multipole moments as 
long as a > 4 x 1CT 4 (with e 1 ~ 10~ 2 ), but starts to 
modify the Ct's at intermediate scales for smaller values 
in accordance with the previous discussion (see FigQJ. 
Let us also recall that a is the dimensionless parameter 
controlling the frequency of the oscillations. It is given by 
the ratio of the Hubble parameter during inflation to the 
scale at which the new physics is supposed to show up, 
<7 = H/M c . Therefore, we conclude that the sampling 
or splinning option has an important effect and must be 
treated with great care. However, despite the above dis- 
cussion, we show below that we essentially recover the 
same fit as the one displayed in Ref. up to some small 
differences. 



B. Lensing 

It has been suggested that CMB lensing effects could 
blur superimposed oscillations. This could have an in- 
fluence on the fit found in Ref. Q] since this one cor- 
responds to high frequency wiggles. However, we have 
checked by mean of the full sky lensing routines imple- 
mented in CAMB jjj that, even in this case, the wiggles 
around the first Doppler peak remain almost unchanged. 
This result is expected since the damping due to lensing 
is mainly proportional to the amplitude of the superim- 
posed oscillations 0- Moreover, one has to keep in mind 
that oscillations do not correspond to a localized feature 
in the Fourier space as studied in Ref. [6(, but, on the 
contrary, constitute a modification which, in some sense, 
is spread everywhere. 



C. Exploring the "fast" parameter space 
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FIG. 1: The TT, TE and EE angular power spectra for the 
low frequency (LF) and high frequency (HF) models. A zoom 
of the Cj T, s in the first Doppler peak region is also shown 
(black curve) and compared with the standard slow-roll pre- 
diction (blue curve) calculated the with same cosmological 
parameters. 



We now turn to the effects that could influence the 
search of the likelihood maxima and/or the analysis of 
their statistical significance. 

As noticed in Ref. the fact that it is necessary 
to boost the computation accuracy in order to correctly 
transfer the oscillations significantly increases the com- 
putation time for one model. The use of the sampling 
(see the discussion above), as in Ref. pj, permits us 
to explore the full (9-dimensional) parameter space us- 
ing Monte Carlo methods implemented in the C0SM0MC 
code 0- Nevertheless, this exploration remains clearly 
limited and may have very well missed the global max- 
imum of the likelihood function. As a consequence, the 
fit found in Ref. [I( should rather been called "a better 
fit" than "the best fit" . 

In the approach advocated here, where the sampling 
and splinning have been removed, the situation becomes 
even worst. In this case, an exploration of the full param- 
eter space, even limited, is no longer possible. In order to 



tackle this problem, we have fixed the cosmological pa- 
rameters to their "standard values" , determined from the 
best fit obtained with a vanilla slow-roll power spectrum. 
Then, the parameter space to be explored becomes much 
smaller and now consists in the so-called "fast parame- 
ters" only, i.e. the slow-roll parameters, the frequency 
and the amplitude of the oscillations. This method al- 
lows to compute the CMB transfer functions only once 
and, as a consequence, decreases significantly the com- 
putation time of each Markov chain. But, clearly, one 
should keep in mind that, if the actual best fit is some- 
where else in the parameter space, this method will miss 
it. 



D. Comparing the fits 

In Ref. 0, the "best fit" obtained after a limited ex- 
ploration of the full parameter space and characterized 
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FIG. 2: Cumulative residual y 2 f° r the LF and HF mod- 
els, compared to the standard slow-roll one. The dashed 
red curve represents the temperature contribution, the dotted 
blue curve the polarization contribution and the black curve 
is the sum of these two effects. 

by x 2 — 1415.4 (1340 DOF) has been compared to the 
best inflationary fit published in Ref . |4| , the % 2 of which 
is given by \ 2 — 1431 (1342 DOF). Hence the number 
Ax 2 — 15 reported in that article. However, the best fit 
of Ref. 4] has been derived under the assumption that the 
primordial power spectrum is of the form fc™s _1 which is 
not exactly the slow-roll prediction In addition, no 
gravitational waves have been included whereas this is 
automatic in the slow-roll approach because of the consis- 
tency check of inflation. Therefore, the number A\ 2 — 15 
does not describe the effect of the wiggles only. In order 
to disentangle the influence of the wiggles from the effects 
of the other parameters, one should first determine the 
best fit obtained with a standard slow-roll power spec- 
trum and then compare this fit to the one obtained after 
the addition of the oscillations, keeping the same values 
for the cosmological parameters. Our best slow-roll fit 
(without the oscillations) corresponds to a model with 
X 2 ~ 1429.3 (1342 DOF), i.e. a difference of A X 2 ~ 2 in 
comparison with the best fit of Ref. |4( . Therefore, the 
effect of adding oscillations is more accurately described 
by Ax 2 — 13 rather than 15. In the following, the evalu- 
ation of Ax 2 will always refer to a comparison with this 
slow- roll fit. 



E. The new improved fits 

In this section, we present two interesting models cor- 
responding to two local maxima of the likelihood, the 
improvements discussed before having been taken into 
account. The results are summarized in Table [I] One 
model, labeled "HF" (high frequency), corresponds to 
very rapid oscillations, while the other, labeled "LF" (low 
frequency), has a smaller frequency. The HF model is es- 
sentially similar to the "best fit" model found in Ref. [lj 
(see the Table I in that article), hence justifying the claim 
that one recovers almost the same solution despite the 
new effects studied in this addendum. The improvement 
of the x 2 is now A^ 2 — 11 with two additional parame- 
ters, that is to say slightly less good (Ax 2 — 2) but still of 
the same order of magnitude as the fit found in Ref. . 
The corresponding F probability is -Fproba — 0.6% and 



hence the improvement remains statistically significant. 
As discussed in Ref. lj, this model suffers from a severe 
back-reaction problem. 

As already mentioned, the LF model seems to be also 
favored by the data. Admittedly, the improvement of the 
X 2 is less good with Ax 2 — 8.5 for two additional param- 
eters but corresponds to -Fproba — 1-87%, demonstrating 
that it is also statistically significant. Furthermore, this 
model no longer suffers from the back reaction problem 
which is, from the theoretical point of view, an impor- 
tant advantage. Indeed, as can be seen from Eq. (11) 
in Ref. pj, back reaction effects are important for small 
values of cr , i.e. in the case of high frequency wiggles. 
This explains why the LF fit can satisfy the back reaction 
constraint. 

Finally, in Fig. ^ we present the HF and LF fits and 
in Fig. |21 we plot the cumulative residual x 2 with respect 
to the standard slow-roll model. 



III. OTHER COSMOLOGICAL OBSERVABLES 

In Ref. 0, we have only presented predictions for the 
quantities CJ T and CJ E . However, two other important 
observables, for which data already exist or will be avail- 
able very soon, are the matter power spectrum and the 
EE polarization channel, Cf E . In this addendum, we 
calculate these two observables in the case where super- 
imposed oscillations are present. 

In Fig. |31 we present the matter power spectrum cor- 
responding to the two HF and LF best fit values given 
in Table together with the current SDSS 01 decon- 
volved power spectra in the linear regime. It is obvious 
that the oscillations are transfered from the initial power 
spectrum to the matter power spectrum since there are 
linked by a transfer function only, P(k) ~ T(k)PQ. For 
the HF and LF fits, we see that the oscillations are well 
within the error bars. We have also computed the convo- 
lution of the LF and HF oscillatory power spectra with 
the sloan survey windows functions (the blue curves in 
Fig-EJ- It is clear that the HF and LF convolved matter 
power spectra are fully degenerate with the vanilla slow- 
roll power spectrum. As a consequence, we conclude that, 
with the currently available large scale structure data, no 
constraint can be put on the parameters controlling the 
shape of the oscillations. 

Finally, in Fig.^ we present the EE polarization mul- 
tipole moments. As can be seen from the figure, the os- 
cillations are also transfered to Cf E . It is likely that the 
future polarization measurements will play an important 
role in deciding whether the superimposed oscillations 
are really present in the data since the drop in the x 2 is, 
for the moment, almost insensitive to polarization (see 
Fig. [21). With more accurate polarization data at our 
disposal, this situation will certainly change. 
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Type 


|at|ff 
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Q, h h 2 






T 


-^scalar 


6l 


e 2 


n s = 1 — 2e 1 — e 2 


X 2 /DOF 


HF 
LF 


0.107 
0.047 


1.7 x 10" 4 
3.4 x 10" 4 


0.734 
0.734 


0.024 
0.024 


0.116 

0.116 


0.74 
0.74 


0.1294 
0.1294 


23.6 x 10" 10 

23.7 x 10" 10 


0.012 
0.010 


-0.0298 
-0.0226 


1.005 
1.002 


1418.4/1340 
1420.8/1340 



TABLE I: Best HF and LF fit parameters from the WMAP data. The pivot scale and the time rj (see Ref. p|) are still chosen 
such that k*/a = M c = 0.01 Mpc" 1 . 
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FIG. 3: Predicted matter power spectra for the LF and HF 
models (black curve), compared to the deconvolved SDSS 
data [To|. The convolved power spectra by the sloan survey 
windows function (blue curve) are also plotted. 



IV. CONCLUSIONS 



the EE polarization multipole moments. With the avail- 
able data, we have concluded that these two observables 
do not affect the results reached in Ref. [jj. 

A last remark is in order here. In Ref. the statistical 
significance of the wiggles has been discussed by means of 
the F test. This was motivated by the fact that the F test 
only requires the knowledge of the likelihood maxima for 
the different models under consideration. However, by 
only varying the "fast parameters", one may expect to 
efficiently probe the corresponding subspace. Therefore, 
an interesting improvement of the present work would 
be to compute the statistical evidence in this subspace, 
thus providing us with a different test of the statistical 
significance of the wiggles. 



In this addendum, we have investigated various ef- 
fects that could influence the presence of wiggles in the 
CMB multipole moments. We have shown that the splin- 
ning is an important parameter which, in the presence of 
high frequency superimposed oscillations, can influence 
the shape of the of CMB angular spectrum at intermedi- 
ate scales. We have shown that the model presented in 
Ref. Q] remains favored by the data (up to a variation of 
Ax 2 — 2). Two other cosmological observables have also 
been presented, namely the matter power spectrum and 
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